Suppression and emergence of granular segregation under cyclic shear 
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Mixtures of granular materials tend to segregate under various disturbances, such as shear and 
gravity. While convective flows have been implicated in many segregation processes, the particle- 
scale rearrangements associated with segregation are not well understood. We have performed 
experiments on a three-dimensional bidisperse mixture in a split-bottom geometry, under both 
steady and cyclic shear. The pile segregates under steady shear, but the cyclically driven system 
either remains mixed or segregates slowly, depending on shear amplitude. The motion of individual 
grains, imaged in three dimensions using a Refractive Index Matched Scanning technique, shows no 
signs of particle-scale segregation dynamics that precede bulk segregation. Instead we find that the 
transition from non-segregating to segregating flow is accompanied by a transition to significantly 
less reversible particle trajectories (even for monodisperse systems), and the emergence of a collective 
convective flow field. 

PACS numbers: 45.70.Mg, 45.70.Qj, 47.57.Gc 
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When mixtures of granular materials are continuously 
disturbed by external forcing, such as vibration, grav- 
ity, rotation, or shear, the grains often separate based 
on their species properties, such as size or density [1- 
[iq| . This phenomenon, known as segregation^ has been a 
subject of scientific interest for several decades because 
of its widespread applicability in both nature and indus- 
try, from stratification of avalanche deposits to the non- 
uniform settling of mixed nuts and cereals. 

Previous studies look at the onset and patterns of seg- 
regation, and several theoretical models have been pro- 
posed [H, '3, |l2[ . A universal feature of the models 
is gradient-driven segregation flux. In the specific case 
of dense shear-driven segregation, however, the relative 
contributions of gradients in friction, shear, and kinetic 
temperature are not well established. In addition, most 
candidate models for shear-driven segregation are linear 
12 . [l3| . A key nonlinearity missing from the current 
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models is granular convection, which can play a vital role 
in the segregation process |3|,|9|,ll^. Indeed, the behavior 
of experimental shear-driven systems deviates from what 
current models predict [iq| . 

In particular, convective flows can lead to the rise of 
large particles in vibrated beds Even under quasi- 
static vibrations with very low kinetic temperature and 
convection, critical size ratios for the rising of large in- 
truders have been observed [l, [HI. While this has in- 
spired further work on vibro-excited systems, there is lit- 
tle discussion in the current literature about shear-driven 
segregation in the absence of kinetic temperature. 

One distinguishing feature of shear-driven flows is that, 
unlike vibro-excited systems, they can exhibit flow driven 
by bulk shear forces, with particle-particle interaction po- 
tential energies that are orders of magnitude larger than 
kinetic energies Driven slowly, granular shear flows 
then can exhibit close to reversible particle trajectories 



of segregation are based on gradients, therefore invari- 
ant under reversal of the shear direction, and predict 
the same segregation fiux for both steady and oscillatory 
shear. Indeed, the absolute size of individual particles 
relative to the amount of shear does not enter into these 
models, nor does a characteristic length scale over which 
gradient-induced segregation manifests. 

In this Letter, we provide insight to these challenges 
by presenting experimental results from a dense three- 
dimensional (3D) shear-driven bidisperse granular ma- 
terial. Using the Refractive Index Matched Scanning 
(RIMS) imaging technique [l9|, we are able to experimen- 
tally determine the trajectories of almost all particles in 
the mixture as the system is slowly sheared, both steadily 
and cyclically, in a split-bottom geometry. We directly 
measure the bulk segregation patterns under steady and 
cyclic shear, and characterize segregation at the micro- 
and mesoscale by quantifying the reversibility of particle 
rearrangements and examining the resulting secondary 
flows of the system. 

Setup and procedure - The granular system is a dense 
bidisperse mixture of spherical Polymethyl methacry- 
late (PMMA) spheres, with small and large diameters 
of Ds = 1/8 inch and Dl = 3/16 inch, respectively. The 
mass ratio between the two species of grains = ^) 
is selected such that the number ratio of small to large 
grains ^ = 2. The sides of the square tank containing 
the pile are 15 cm long and the pile height is 4.6 cm. 
The pile height is approximately the radius of the split- 
bottom shearing disk, Rs = 4.5 cm, as well as IOI^l- 

The grains are immersed in an index-matched inter- 
stitial fluid, Triton X-100. The index-matching allows 
light from two laser sheets, mounted on either side of the 
tank, to pass through undefiected. We also add a laser 
dye, Nile Blue 690 Perchlorate, that fluoresces at the laser 




FIG. 2. Cross-sections of the granular pile above the shearing 
disk (a) in initial mixed state (b) after 29.3 revolutions of the 
shearing disk. These particular cross-sections are taken at 3.8 
cm from the bottom, or, about 2 large grain diameters from 
the free surface of the pile. Red and blue circles represent 
large and small grains, respectively. 




FIG. 1. (a) The split-bottom shear cell used in this experi- 
ment, with a sample horizontal cross- section (the second laser 
sheet on the opposite side of the tank and a high sensitive 
camera that captures images from the top down are not pic- 
tured), (b) Sample trajectories of a small (green to blue) and 
large (yellow to red) grain under steady shear. The black 
circle represents the shearing disk. 



wavelength (635 nm), as well as a small amount of hy- 
drochloric acid to stabilize the mixture. With the laser 
sheets aligned at a particular height, a high-sensitivity 
Sensicam camera captures a horizontal cross-section of 
the pile from the top-down. The laser sheets and cam- 
era move incrementally along stepper motors, with cross- 
sections taken about every 200 /im in height. All of the 
cross-sections together construct a full 3D image of the 
pile. An illustration of our setup is shown in Figure [TJa). 

We consider the response of the granular system to 
both steady and cyclic shear, with amplitudes of 10 and 
40 degrees. The shearing disk rotates at a rate of 1 
mrad/s, with static 3D images captured every 2 degrees. 
At this shearing rate, the flow profiles of the pile resemble 
that of a dry pile and are also rate- independent [20]. 

When an experiment is complete, particle center po- 
sitions are determined using a 3D-adapted convolution 
kernel [2l|, which also distinguishes between large and 
small grains. In each frame, we are capable of extracting 
the position of at least 95% of particles, with a particle 
center resolution of about 100 /im. Then, a Lagrangian 
predictive particle tracking algorithm [22| identifies indi- 
vidual particles and their trajectories. Examples of in- 
dividual particle trajectories, for a large and small grain 
under steady shear, are shown in Figure [TJb). 

Bulk Segregation - A standard segregation pattern that 
is seen in many polydisperse mixtures is the Brazil Nut 
Effect (BNE), in which large particles rise to the top of 



the pile, leaving smaller sizes to settle toward the bot- 
tom. In this slowly sheared system, a pattern similar to 
the BNE is observed under steady shear. Figure [2] shows 
how a cross-section at a height of 3.8 cm (8 large grain 
diameters) becomes primarily populated with large par- 
ticles after almost 30 revolutions of steady shear. This 
height is close to the free surface of the pile, only two 
large grain diameters away. This effect is also indicated 
by the sample trajectories shown in Figure [TJb). The pile 
as a whole rearranges such that large particles inhabit a 
cylindrical region at the top 60% (~ VL/Vtotah where Vl 
is the total volume of large grains and Vtotai is the to- 
tal volume of all grains) of the pile, directly above the 
shearing disk. To quantify bulk segregation, we calculate 
the volume fraction of large and small grains in this re- 
gion of interest. Figure [3{a) illustrates that over about 
30 rotations of the shearing disk, the pile is continually 
segregating. 

Under oscillatory shear, however, segregation is not 
necessarily observed. The results for bulk segregation 
under steady shear and the two oscillatory amplitudes 
considered are summarized in Figure [Ifb). When the 
amplitude is 10 degrees, the volume fractions of large and 
small particles appear to remain constant throughout. 
This indicates that the system is not segregating over the 
time interval measured. Individual particle trajectories 
are consistent with this observation, as discussed in the 
following section. 

However, for sufficiently large amplitudes of oscillatory 
shear, the pile does slowly segregate. For an amplitude 
of 40 degrees, the rate of segregation is very small, but 
noticeable over the measured time period. As seen in Fig- 
ure [3jb), there is a clear distinction in the evolution of 
a segregating system (steady shear and 40 degree cyclic 
shear) and a non-segregating oscillatory system (10 de- 
gree cyclic shear). 

Microscale reversibility - The RIMS technique allows 
us to also probe segregation microscopically by consider- 
ing particle contact networks and mean square displace- 
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FIG. 3. (a) Volume fractions of large, small, and all grains 
in the top 60% of the pile, right above the shearing disk, 
(b) Same as (a), but with large and small grain volume frac- 
tions under all shearing procedures considered: steady (red 
and blue). Or = 10° (orange and cyan), and Or = 40° (gold 
and green). Inset: the average rate of change of of large 
and small grains under the considered shearing procedures 
for |A6'| > 1000°. 



merits (MSDs) within the shear zone of the split-bottom 
geometry, where most rearrangements occur [lH, [23, 24 1. 

In order to verify that the granular material is segre- 
gating, we must see that grain motion is irreversible in 
terms of both grain-to-grain contacts and displacements. 
The pile starts in a highly mixed state, so each individ- 
ual grain touches a combination of small and large grains. 
Under segregation, particle contacts must be constantly 
changing, in order to form new contacts with similar par- 
ticles. At the same time, segregation requires significant 
and irreversible motion within the pile, so that a highly 
stratified state as described by the BNE can be obtained. 

We determine grain contacts by finding particle pairs 
that, within some distance cutoff, are increasingly likely 
to move perpendicular to each other by sliding or rolling, 
as described by Herrera, et al [25]. Three contact length 
cutoffs are defined: large-large, large-small, and small- 
small. MSDs are simply calculated from the average dis- 
placement of grains from an initial reference state to a 
new state at a later time. 

Figure lU^a) shows the fraction of broken links at the 
end of a single shear cycle, with respect to various ref- 
erence networks taken at the end of the previous cycle. 
For all three contact types, there is a decreasing fraction 
of contacts that are broken over time for the 10 degree 
experiment; however, the fraction of broken links for the 
40 degree experiment remains fairly steady from cycle to 
cycle. Similar trends are observed if we measure MSDs 
of large and small grains over individual cycles, as shown 
in Figure m^b). This distinguishing behavior of reversible 
and irreversible flows was also seen by Slotterback, et al 
in a monodisperse system with particle size Dl^ with a 
similar critical shear amplitude y^]. This indicates that 
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FIG. 4. (a) Fraction of broken grain contacts and (b) MSDs of 
large and small grains over a single cycle versus that particular 
cycle number. 



in this split-bottom geometry, there is a connection be- 
tween a shear-driven granular material undergoing irre- 
versible flows and having the capacity to segregate. 

We have described the bulk segregation patterns 
we see under steady and cyclic shear, and linked 
our segregating/non-segregating regimes to changes in 
particle-scale reversibility. To paint a more complete pic- 
ture of how the granular material rearranges, we now 
connect both bulk behavior and particle-scale motion to 
the secondary flow profiles for our bidisperse system. 

Secondary flows - The motion of the grains is primarily 
in the (/)-direction, as grains near the surface are driven 
locally by the rotation of the shearing disk. However, the 
grains are also free in move in the other two cylindrical 
coordinates, r and z. The time- and azimuthal- averaged 
motions of grains in r and z are referred to as secondary 
flows. Simulations in the split-bottom geometry suggest 
that secondary flows are a key component of the segre- 
gation process [9]. 

For steady shear, the secondary flows for large and 
small particles are shown in Figure [5f a). The flow pro- 
files of the two species are vastly different; while the large 
particles form two distinct convection rolls, the small par- 
ticles primarily drift toward the bottom. There is a net 
downward flux of small particles above the shearing disk, 
which is in line with the bulk segregation pattern. Also 
shown in Figure [5] is the outer border of the shear zone, 
which is determined from the 2;-dependence of the az- 
imuthal flow profile [23I, (2^ . The secondary flows exist 
within this boundary, suggesting that the segregating re- 
gion is necessarily determined by the size of the shearing 
disk. 

For oscillatory shear, the secondary flow profiles are 
dependent on the shear amplitude. For an amplitude 
of 10 degrees, the secondary flows of both species are 
suppressed, as illustrated in Figure [5fb). Due to the re- 
versibility of the system, this is to be expected. When the 
shear amplitude is increased to 40 degrees, granular con- 
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FIG. 5. Time- and 0- averaged flows of large and smafl grain flows plotted together for all shearing procedures considered: (a) 
steady, (b) Or = 10°, (c) Or = 40°. As a guide, the outer boundary of the shear zone, determined from the azimuthal flow 
proflles, is also shown. All vectors are on the same arbitrary scale. 



vection remains mostly undeveloped for large particles. 
However, the small particles, particularly those toward 
the bottom of the pile and close to the shearing disk, 
drift downward, as shown in Figure [5fc). For higher am- 
plitudes, as grain irreversibility increases, we expect to 
see stronger downward drifts of small particles through- 
out the pile, as weh as granular convection of large grains. 
For this particular amplitude, we see the lack of convec- 
tion as a contribution to the diminished segregation rate. 

Discussion - Using the RIMS technique, we image and 
track individual grains in a dense 3D bidisperse granular 
mixture under the influence of shear in a split-bottom 
geometry. We observe that in a regime of slow, rate- 
independent steady shear, the system segregates with 
large particles occupying the top portion of the pile di- 
rectly above the shearing disk. When the system is 
sheared cyclically, there is a critical shear amplitude be- 
low which segregation is suppressed. 

We find that bulk segregation may be linked to micro- 
scopic reversibility of small and large grain trajectories. 
Under 10 degree oscillatory shear, small and large grains 
are increasingly reversible over time, as observed struc- 
turally from their grain contacts and spatially from their 
MSDs. However, at 40 degree oscihatory shear, the de- 
gree to which the large and small grains are reversible 
remains steady over all cycles. These trends of micro- 
scopic reversibility are very similar to observations of the 
monodisperse system [18], indicating that they are not 
due to any additional forces, unique to granular mixtures, 
that drive segregation. 

Finally, we investigate the secondary flow profiles 
brought about by the (ir)reversibility of the system. Un- 
der steady shear, we find clear convective flows of large 
grains, which help drive the segregation mechanism by al- 
lowing large grains to settle at the top of the pile, while 
there is a net flux of small grains falling between par- 
ticle gaps toward the bottom of the container. These 
flows are also confined within the outer boundary of the 
shear zone, suggesting that the size and shape of these 
fiows are dependent on the system geometry. For small 
amplitude oscillatory shear, however, these fiows are sup- 



pressed. Just past the critical shear amplitude for seg- 
regation, convection of large grains has not been allowed 
to completely form, but there is a clear downward drift 
of small grains near the shearing disk. 

In this system, we observe distinct rates of segregation 
for steady and adequately large amplitude oscillatory 
shear. However, we cannot immediately attribute these 
rates to relevant system parameters, such as shear ampli- 
tude, diffusion coefficients, or gradients of local strain. A 
recent study also indicates the presence of gravity is a de- 
terminant of granular convection [26] , which would be rel- 
evant for segregation rates in astrophysical applications. 
Future work on this system should focus on probing the 
functional dependence of segregation rate on these vari- 
ables, which is key to determining a general predictive 
model for shear-driven segregation [27|. 

We are encouraged by the observation of a critical 
shear amplitude that brings about segregation. In fact, a 
critical strain of similar magnitude also exists for fracture 
[25] and irreversibility [18], suggesting that this granular 
system has a universal length scale that drives inherently 
bulk rearrangements, yet arises directly from particle- 
scale dynamics. This length scale could be instrumental 
in the development of not only segregation models, but 
also predictive constitutive models for granular fiows that 
rely on a characteristic length for cooperative rearrange- 
ments [28]. 

Acknowledgments - We would like to acknowledge 
enlightening discussions with Mitch Mailman, Kerstin 
Nordstrom, and Anthony Thornton, as well as technical 
and analytic assistance from Steven Slotterback. Finan- 
cial support for this work came from the National Science 
Foundation (DMR0907146) and the Defense Threat Re- 
duction Agency (HDTRAl-10-0021). 



mjharrin@um d.edu| 
^ I wlosert @um d . edu 
[1] S. B. Savage and C. K. K. Lun, 



5 



|J. Fl uid Mech. 189, 311 (1988) 
J. Duran, J. Rajchen bach, 



and E. Clement, 



[Phys. Rev. Lett. 70, 2431 (1993)| 
J. B. Knight, H. M. Jaeger, 
[Phys. Rev. Lett^ 70, 3728 (1993)1 
D. Khakhar, J. McCarthy, 
Phys. Fluids 9, 3600 (1997) 

J_. Gray and 

"Cont . Mech. Therm. 9, 341 (1997) | 
N. J ain, J. M. Ottino, a nd 
Phys. Rev. E 71, 051301 (2005) 
C. R. Charles, Z. S. Khan, 
Gran . Matt. 8 , 1 (2006) 
A. Golick 



and S. R. Nagel, 



and J. Ottino, 



K. 



Hutter, 



R. M. Lueptow, 



and S. W. Morris, 



L. 



and 



K. 



Phys. Rev. E 80, 042301 (2009)1 



E. Daniels, 



Y. Fan and K. M. Hih, Phy s. Rev. E 81, 041 303 (2"oToy. 
L. B. H. May, L. A. Golick, K. C. Phillips, M. Shearer, 
and K. E. Daniels, Phys. Rev. E 81, 051301 (2010), 
V. Garzo and J. W. Dufty, Phys. Fluids 14, 1476 (2002) , 



Gray and A. 

|Proc. R. Soc. A 461, 1447 (2005)[ 
Y. Fan and K. 



M. 



New Journal of Physics 13, 095009 (2011), 
R. Khosropour, J. Zirinsky, H. K. Pak, 
Behringer, Phys. Rev. E 56, 4467 (1997), 
S. Dippel and S. 

[J. Phys. I France 5, 1527 (1995)] 
L. Kondic and R. P. 



Thornton, 
Hih, 
and R. P. 

Luding, 
Behringer, 



[17; 
[is: 

[19 

[20; 

[21 

[22; 

[23' 
[24 

[25; 

[26; 

[27 
[28; 



[Euro. Phys. Lett. 67, 205 (2004)', 
D. Pin e, J. Go llub, J . Brady, 



[Nature 438, 997 (2005) 



S. 
M. 



Slotterback, 
van Hecke, 



M. 
M. 



Mailman, 
Girvan, 



and A. Leshansky, 

K. Ronaszegi, 
and W. Losert, 



Phys. Rev. E 85, 021309 (2012), 

J. A. Dijksman, F. Rietz, K. A. Lorincz, M. van Hecke, 
and W. Losert, [Rev. Sci. Inst r. 83, 011 301 (2012), 
J. A. Dijksman, E. Wandersman, S. Slotterback, C. R. 
Berardi, W. D. Updegraff, M. van Hecke, and W. Losert, 
[Phys. Rev. E 82, 060301 (2010) [ 

J.-C. Tsai and J. P. Gollub, 

[Phys. R ev. E 70, 031303 (2004) , 

N. T. Ouellette, H. Xu, and E. Bodenschatz, 
|Exp. Fluids 40, 301 (2006) , 

D. Fenistein and M. van Hecke, "Nature 425 , 256 (2003) 

D. Fenistein, J. W. van de Meent, and M. van Hecke, 
|Phys. R ev. Lett. 92, 094301 (2004) , 

M. Herrera, S. McCarthy, S. Slotterback, 

E. Cephas, W. Losert, and M. Girvan, 
[Phys. Rev. E 83, 061303 (2011)[ 

N. Murdoch, B. Rozitis, K. Nordstrom, S. F. Green, 

P. Michel, T.-L. de Lophem, and W. Losert, 

[ Phys. Rev. Lett . 110, 018307 (2013), 

A. Thornton, T. Weinhart, S. Luding, and O. Bokhove, 

Jnt. J. Mod. Phy s. C 23, 1240 01^(2012) , 

K. Kamrin and G. Koval, 

[Phys. Rev. Lett. 108, 178301 (2012)] 



